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Abstract. We have carried out a sensitive search for 
young pulsars associated with supernova remnants using 
the 76-m Lovell radio telescope at Jodrell Bank. The ob- 
servations were made at 606 MHz using a system with 
a bandwidth of 8 MHz and a set noise temperature on 
cold sky of about 50 K. The survey targeted 33 rem- 
nants in the northern hemisphere and achieved a nominal 
sensitivity of ~ 1 mJy in most cases. Two pulsars were 
discovered in the course of this survey and the known 
pulsar PSR B1952+29 was detected. The new pulsars, 
J0215+6218 and J1957+2831, were found during searches 
of the supernova remnants G132.7+1.3 and G65. 1+0.6 re- 
spectively. Based on a statistical analysis of the present 
sample of proposed pulsar-supernova remnant pairs, we 
conclude that at most 17 associations are likely to be real. 
We find no strong evidence for a genuine association be- 
tween either of the two newly discovered pulsars and their 
target supernova remnants. 

Key words: pulsars: individual (PSRs: 

J0215+6218, J1957+2831) — supernova remnants: indi- 
vidual G132. 7+1.3, G65. 1+0.6 



1. Introduction 

In the standard model for pulsar evolution (Gunn & Os- 
triker 1970), young rapidly rotating pulsars are expected 
to be harboured within the remnants of the supernova ex- 
plosions in which they were formed. The two best-known 
examples of young pulsars, viz. Crab and Vela, were dis- 
covered in the early days of pulsar astronomy (Staelin 
& Reifenstein 1968; Large et al. 1968). The overwhelm- 
ing evidence in favour of their association with the Crab 
and Vela supernova remnants helped to establish the ba- 
sic rotating neutron star model to explain the pulsar phe- 
nomenon (Pacini 1968; Gold 1968). In this paper, we de- 
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scribe a survey carried out to find pulsars associated with 
supernova remnants in the northern hemisphere. 

In the past few years, the number of proposed as- 
sociations between pulsars and supernova remnants has 
risen to around 30 (see Kaspi 1996 for a recent review). 
This is a result of a number of different approaches: 
high frequency searches of the Galactic plane (Clifton et 
al. 1992; Johnston et al. 1992); cross-correlations of the 
pulsar and supernova remnant catalogues (Caraveo 1993; 
Fiirst et al. 1993; Kulkarni et al. 1993); targeted searches 
for pulsars in remnants (Frail & Moffett 1993; Gorham et 
al. 1996; Kaspi et al. 1996) and remnants around young 
pulsars (Frail et al. 1994). 

The number of spurious associations in this sample, i.e. 
those cases where the pulsar and the supernova remnant 
are merely in chance alignment on the sky, is likely to be 
significant. Theoretically, the expected number of associ- 
ations has a complex dependence on a number of factors: 
the luminosity function and initial spin period of young 
pulsars (Narayan 1987; Narayan & Schaudt 1988), their 
beaming fraction (Lyne & Manchester 1988; Biggs 1990), 
space velocities (Shull et al. 1989; Caraveo 1993; Lyne 
& Lorimer 1994), the evolution of supernova remnants 
(Gaensler & Johnston 1995a&b) as well as the fraction 
of all supernovae that produce pulsars. In a detailed sta- 
tistical analysis, Gaensler & Johnston (1995b) concluded 
that fewer than 10 associations are expected to be real. 

Statistical studies require accurate flux density upper 
limits for any pulsars undetected in remnants. Our search 
complements recent targeted searches on southern super- 
nova remnants at Parkes (Kaspi et al. 1996), on those rem- 
nants visible from Arecibo (Gorham et al. 1996), as well 
as an earlier survey at Jodrell Bank (Biggs & Lyne 1996). 
The present survey was successful, discovering two new 
pulsars as well as substantially improving on previous up- 
per limits for the flux density of any pulsars in the target 
remnants. 

The plan for this paper is as follows: In Sect. 2 we 
describe the survey observations and data reduction tech- 
niques. In Sect. 3 we estimate the sensitivity of the sur- 
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Table 1. 33 supernova remnants targeted by the survey. From left to right the columns give the remnant name based on its 
Galactic coordinates, any alias(es) by which the remnant may be called, the equatorial coordinates, the approximate angular size, 
a classification tag (S - shell remnant, ? - unknown type, F - filled centre). Question marks denote poorly- known observational 
quantities. All these data are taken from Green's (1996) catalogue of supernova remnants to which the interested reader is 
referred for further details. We also list the estimated contribution to the system temperature from the sky background (TLky) 
and the remnant itself (T TCU1 ), the number of telescope pointings used to observe each remnant (-/Vp nt ), as well as the estimated 
minimum flux density at 606 MHz (S m i n ), above which a pulsar would have been detected by the survey (see text). 



vey. We summarise the results of the survey in Sect. 4. In 
Sect. 5 we discuss the validity of the association between 
the detected pulsars and their target remnants and com- 
ment on the implications of the flux density upper limits 
for the luminosity of pulsars at birth. In Sect. 6 we carry 
out a simple model-free statistical analysis to deduce the 



number of associations that would be expected by chance. 
Finally in Sect. 7 we present our conclusions. 

2. Survey Observations 

In order to complement recent searches in the southern 
hemisphere and from Arecibo, we chose to search primar- 
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ily those supernova remnants with declinations north of 30 
degrees. Our final sample of 33 supernova remnants was 
selected from Green's (1996) catalogue and is summarised 
in Table 1. 

The survey observations were made with the 76- 
m Lovell radio telescope operated by the University of 
Manchester at Jodrcll Bank, on 5 separate observing ses- 
sions between 1994 June and 1996 April. The centre fre- 
quency for the observations was 606 MHz with a total 
bandwidth of 8 MHz. For those remnants larger than the 
telescope beam width (0.5 degrees FWHM) we covered the 
full area of the remnant using up to 13 individual telescope 
pointings (see Table 1). Each telescope pointing was of 35 
minutes duration. 

The incoming radiation was split by the feed into two 
orthogonal linear polarisations and amplified by a pair 
of cryogenically-cooled field effect transistors. The signals 
then entered wide-band filters and, following a second 
stage of amplification, were mixed down to an intermedi- 
ate frequency centred on 8 MHz before entering a pair of 
64 x 0.125 MHz filtcrbanks. The resulting 64 polarisation 
pairs were then detected, summed, integrated and one- 
bit digitised every millisecond. The data were then passed 
to an on-line VAX computer which wrote the data to a 
standard Exabyte magnetic tape for off-line processing. 

The off-line search for periodic signals in the data con- 
sisted of a Fourier analysis of time sequences of dedis- 
pcrsed data in order to search for significant features in the 
power spectrum. This was followed by a time series analy- 
sis to optimise the parameters of the most promising spec- 
tral features. The complete search procedure is very simi- 
lar to that used for the highly successful Parkes southern 
sky survey described by Manchester et al. (1996). Briefly, 
the data were de-dispersed and Fourier transformed for 90 
trial values of dispersion measure (DM) between zero and 
995 cm~ 3 pc. For any periodic signal with a small duty cy- 
cle (5-10% is typical of most pulsars), the resulting power 
spectrum from the Fourier transform consists of a family 
of harmonic spikes with the fundamental corresponding to 
the signal frequency. To increase the sensitivity to pulsars 
with narrow pulses, higher order harmonics can be added 
onto the fundamental (Lyne 1988). For example, to add 
all 2 nd harmonics onto their corresponding fundamentals, 
we stretch the lower half of the amplitude spectrum by 
a factor of two and add this to the original unstretched 
spectrum. While the effect of this harmonic summing pro- 
cess increases the noise by a factor of y/2, the amplitude 
of the signal may increase by more than a factor of y/2 
giving a net increase in the signal-to-noise ratio. By re- 
peating this process several times we effect a search in 
pulse duty cycle. In our analysis the harmonic summing 
was performed 4 successive times, i.e. producing power 
spectra summed over 1, 2, 4, 8 and 16 harmonics. During 
this stage periodic signals from known sources of terres- 
trial interference, e.g. the 50 Hz mains power line, were 
excised from the data by zeroing the appropriate portions 



of the power spectrum. Due to strong amounts of terres- 
trial interference in the low-frequency part of the power 
spectrum, we did not consider signals with frequencies be- 
low 0.2 Hz. This degraded our sensitivity only slightly to 
pulsars with periods greater than 5 seconds which could 
still be detected in their higher order harmonics. 

For each value of DM, the largest features in the power 
spectra and harmonic sums were then sorted so as to pro- 
duce a list of "pulsar suspects" . All suspects with spectral 
signal-to-noise ratios greater than 7 were then folded in 
the time domain for a range of periods and DMs around 
the nominal values found in the first stage in an attempt 
to optimise these parameters and the signal-to-noise ra- 
tio. The output from this stage was visually inspected 
(see also Fig. 5 of Manchester et al. 1996) and the best 
suspects were stored for rc-obscrvation. 

3. Search Sensitivity 

The sensitivity of the survey can be expressed in terms 
of the minimum detectable flux density S m i n which is a 
function of a number of parameters. Following Dewey et 
al. (1984), we write: 

_ arjT^ ( W \ 1/2 

In this expression a is the threshold signal-to-noise ra- 
tio above which a detection is considered significant (7 in 
our case), r\ is a constant ^1.3 which reflects losses to 
hardware limitations, G is the gain of the telescope (1.0 
K Jy _1 for the Lovell telescope operating at 606 MHz), n 
is the number of polarisations used, Av is the observing 
bandwidth, T sys is the system temperature, t is the inte- 
gration time, P is the period of the pulsar and W is the 
observed width of the pulse. With the parameters for this 
survey, the above expression simplifies to 

/ W \ 1/2 

S min ~0.05T sy J — — J mJy. (2) 

The system temperature T sys is the sum of separate 
components: the set noise of the receiver T sct ; the sky 
background noise T s k y and the contribution T rom from 
continuum flux of the supernova remnant. Regular cali- 
bration measurements made during the survey indicated 
T sct to be typically 50 K. The contribution to T sys from the 
Galactic background and the supernova remnant covers a 
large range. This is shown in Table 1 where we list the 
estimated sky background temperatures from a machine- 
readable version of the Haslam et al. (1982) all-sky sur- 
vey, scaled to 606 MHz assuming a spectral index of -2.7 
(Lawson et al. 1987) together with the expected contribu- 
tion from each supernova remnant. To calculate the lat- 
ter values, we estimated the flux density of the remnant 
at 606 MHz from spectral information in Green's (1996) 
catalogue and multiplied this by the beam filling factor, 
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defined as the lesser of unity and (FWHM/I?snr) 2 , where 
-Dsnr is the angular diameter of the supernova remnant. 
Using these values in equation 2, we have estimated the 
minimum flux density required to detect a 0.1 s pulsar 
with a duty cycle of 4% in each remnant. These limiting 
flux densities are listed in Tabic 1 for reference. With typ- 
ical values of ~ 1 mJy at 606 MHz, they demonstrate the 
excellent sensitivity of the survey. We note from Table 1 
that T rcm is a significant factor for only 3 of the 33 SNRs 
searched. 

The observed pulse width W in Eqs. 1 and 2 is likely 
to be greater than the intrinsic width W- lrA emitted at the 
pulsar because of the scattering and dispersion of pulses 
by free electrons in the interstellar medium, and by the 
post-detection integration performed in the receiver. The 
observed sampled pulse profile will therefore be the con- 
volution of the intrinsic pulse width and broadening func- 
tions due to dispersion, scattering and integration and can 
be estimated approximately from the following quadrature 
sum: 

W 2 = Wf nt + 4 mp + tl M + t s 2 catt , (3) 

where i sam p is the data sampling interval, tuM is the 
dispersion broadening across one filterbank channel and 
iscatt is the interstellar scatter broadening. Pulse scatter- 
ing becomes particularly important when observing dis- 
tant pulsars towards the inner Galaxy at frequencies < 
1 GHz. Many of the supernova remnants in our sample 
are relatively nearby ^ 5 kpc so that we do not expect a 
significant effect on our sensitivity due to scattering. 

The effects of sampling and dispersion do, however, 
significantly affect the search sensitivity at short pulse pe- 
riods. This is shown in Fig. 1, where S m i n is plotted against 
P for a pulsar with a dispersion measure of 150 cm~ 3 pc, 
a typical value for a pulsar at the distance of one of the 
remnants in our sample (see Sect. 5.2). We have also cal- 
culated the reduction in sensitivity at short periods due 
to the loss of higher order harmonics in the power spec- 
trum - shown by the abrupt jumps in the sensitivity curve 
shown in Fig. 1. 

It is worth noting that our search had only limited sen- 
sitivity to any high- velocity pulsars that may have moved 
outside the projected boundaries of their parent supernova 
remnants. The dotted curve in Fig. 1 gives the approxi- 
mate sensitivity to a pulsar lying 0.25 degrees outside the 
remnant boundary and indicates a limiting 600 MHz flux 
density ^ 20 mJy. Such bright pulsars should have been 
found in previous large scale surveys of the northern sky. 
Future surveys targeted specifically outside the periphery 
of the remnants will substantially improve the sensitivity 
to fainter pulsars. 




Period (seconds) 

Fig. 1. The limiting sensitivity of the survey as a function 
of pulse period for a pulsar with a DM of 150 cm~ 3 pc and a 
duty cycle of 4% is shown by the solid curve. The calculation 
assumes a total system temperature of ~ 95 K, a typical value 
for supernova remnants observed during the survey. The dotted 
curve shows the approximate sensitivity to a pulsar lying 0.25 
degrees outside the boundary of the remnant (see text). 



4. Results and follow-up observations 

Twenty-five of the most promising pulsar candidates were 
re-observed in April 1996 and the detection of three pul- 
sars was confirmed, B1952+29 and two that were previ- 
ously unknown: PSRs J0215+6218 and J1957+2831. No 
further previously known pulsars were expected to be 
above our detection threshold. The new pulsars lie within 
the boundaries but towards the outer regions of the tar- 
get supernova remnants G132.7+1.3 and G65. 1+0.6 re- 
spectively. No further pulsars were discovered during these 
observations. 

In order to accurately determine the astrometric and 
spin parameters of the two newly discovered pulsars, reg- 
ular follow-up observations at 606 and 1400 MHz have 
been made ever since their discovery as part of the Jodrell 
Bank pulsar timing programme (see for example Shemar 
& Lyne 1996). We have been able to model the pulse ar- 
rival times for both pulsars using standard pulsar timing 
techniques (Manchester & Taylor 1977). Timing solutions 
obtained from these data are given in Table 2. 

Total intensity pulse profiles of the newly discovered 
pulsars at 606 and 1412/1418 MHz are shown in Fig. 2. 
Mean flux densities at both frequencies were determined 
in an identical fashion to that described by Lorimer et 
al. (1995). We also list the characteristic age r c = P/2P 
and the surface magnetic field B 3.2 x 10 19 (FP) 1 / 2 
Gauss which assume a constant dipolar magnetic field and 
short initial period (Manchester & Taylor 1977). 
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Table 2. Observed and derived parameters of the two newly 
discovered pulsars. Quoted errors refer to the last significant 
digit. 




Fig. 2. Total intensity integrated pulse profiles for PSRs 
J0215+6218 and J1957+2831. In each case the whole of 
the pulsar period is displayed. These profiles are freely 
available at the on-line database of pulse profiles 
http:/ /www. mpifr-bonn.mpg.de/pulsar/data. 

5. Discussion 

5.1. Are the detected pulsars really associated with their 
target supernova remnants? 

The previously known pulsar PSR B1952+29 was detected 
during observations of G65. 1+0.6. The pulsar lies outside 
the boundary of the remnant and is clearly not associated 
with G65. 1+0.6 because of its extremely large character- 
istic age (4 x 10 9 yr) and the proper motion measurement 
(Lyne et al. 1982), which shows that it is moving towards 
rather than away from the remnant centre as is required 
for a genuine association. 

The evidence that the newly discovered pulsars are 
associated with their target remnants is not clear. PSR 
J0215+6218 has a period of 549 ms and a DM of 84 cm~ 3 
pc. The distance to this pulsar estimated from its DM and 



Galactic coordinates using the Taylor & Cordes (1993) 
electron density model is 3.2 ± 1.0 kpc, consistent with 
the distance to G132. 7+1.3 of 2.2 + 0.2 kpc (Routledge et 
al. 1991). Timing measurements show that the character- 
istic age of PSR J0215+6218 is 13 Myr, anomalously large 
by comparison with the age of G132.7+1.3, estimated to 
be between 30,000 and 50,000 yr (Leahy et al. 1985). 

In the case of the 308 ms pulsar J1957+2831, there 
is unfortunately no independent distance estimate to the 
target remnant, G65. 1+0.6. However the angular size of 
the remnant (~ 1 degree) indicates that it is likely to be 
closer than the distance of 7.0 ± 2.3 kpc inferred from the 
dispersion measure of PSR J1957+2831 (139 cm" 3 pc). 
Whilst the age of G65. 1+0.6 is also not known, it is likely 
to be much smaller than the characteristic age of PSR 
J1957+2831 (1.6 Myr). 

Thus, in both cases, the pulsar characteristic age is 
anomalously large by comparison with the expected ages 
of the supernova remnants. This suggests that they are 
either not associated with the remnant, or the character- 
istic ages are anomalously large, due perhaps to the initial 
spin period of the pulsar being similar to its presently ob- 
served value. We note in passing that the spectral index 
of PSRs J0215+6218 and J1957+2831 are -1.2 and -0.9 
respectively, typical of many other young pulsars (Lorimer 
et al. 1995) and suggesting that these pulsars may indeed 
be younger than their characteristic ages would suggest. 
Although this remains a possibility, as we shall show in 
Sect. 6, on statistical grounds, neither of the newly dis- 
covered pulsars are likely to be associated with the target 
supernova remnants. 

5.2. The luminosity of pulsars at birth 

In a study of pulsar population statistics Lorimer et 
al. (1993) suggested that there may be no need for a signif- 
icant number of pulsars to be born with 400 MHz radio lu- 
minosities below 30 mjy kpc 2 . Deep surveys of supernova 
remnants, like the present one, can in principle be used to 
test this hypothesis by combining the flux density upper 
limit (S'min) given in Table 1 with the distance to each rem- 
nant (dsNR.) to estimate the minimum luminosity (L mm ) 
that a young pulsar would need in order to be detectable. 
In Table 3 we list the 17 supernova remnants in our sam- 
ple which have reliable distance estimates. Note that we 
have chosen not to use the surface brightness-angular size 
(E — D) relationship (Clark & Caswell 1976) as a means 
of estimating (isNR since it has subsequently been shown 
to be unreliable (Berkhuijscn 1987; Green 1991). 

For each supernova remnant in Table 3, we list the min- 
imum detectable luminosity L min — 2S m i n <i 2 ;NR , where the 
factor of 2 in this expression scales S'min defined for this 
survey at 606 MHz to 400 MHz assuming a typical pulsar 
spectral index of -1.6 (Lorimer et al. 1995) . We find the 
median value of L m i n to be 22 mjy kpc 2 , with 12 of these 
values lying between 1 and 30 mjy kpc 2 . In this sample 
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Table 3. The supernova remnants targeted by the survey with 
published distance estimates. From left to right the columns 
give the remnant name based on its Galactic coordinates, 
alias(es) by which the remnant may be called, the distance 
and reference tag, and the corresponding minimum 400 MHz 
luminosity for a pulsar to be detectable in our survey (see 
text). The references are: a. Lozinskaya et al. (1993); b. Green 
& Gull (1989); c. Landecker et al. (1980); d. Green (1989); 
e. Feldt & Green (1993); f. Tatematsu et al. (1990); g. van 
den Bergh (1971); h. Reich & Braunsfurth (1981); i. Hailey & 
Craig (1994); j. Pineault et al. (1993); k. Albinson et al. (1986); 
1. Green (1996); m. Roberts et al. (1993); n. Routledge et 
al. (1991); o. Reich et al. (1992); p. Landecker et al. (1989); 
q. Routledge et al. (1986) 



of 12 supernova remnants with well defined distances, we 
detected one pulsar (PSR J0215+6218) which has a 400 
MHz luminosity ~ 80 mJy kpc 2 . If we assume that each 
of these remnants contains a radio pulsar born with a lu- 
minosity above 30 mJy kpc 2 , then we expect to see 12/ 
pulsars, where / is the mean beaming fraction. A con- 
sensus on the beaming fraction of young pulsars has yet 
to be established. Frail & Moffet (1993) obtained a value 
/ <~ 0.6 based on a deep imaging search at the VLA. In 
this case we would expect to see ~ 7 pulsars. On the other 
hand, Tauris & Manchester (1997) recently claim / to be 
as low as ~ 0.1 implying a detection of only one pulsar. 
Thus, although our results are consistent with few pulsars 
being born with luminosities below 30 mJy kpc 2 , they are 
still hampered by small number statistics and we conclude 
that yet deeper searches are required to provide a larger 
sample with which to test this hypothesis (see also Kaspi 
et al. 1996). 



5.3. Radio-quiet neutron stars 

The above discussion does not, of course, exclude the pos- 
sibility that neutron stars exist in the remnants which 
do not produce significant radio emission. One classic 
example is the 6.9 s X-ray pulsar IE 2259+586 which 
has long been postulated to be associated with CTB 109 
(Fahlman & Gregory 1981). Despite extensive efforts, no 
radio counterpart to the pulsar has ever been detected. 
Coe et al. (1994) used the VLA at 1489 MHz to set a 3 a 
upper limit to continuum emission of 50 /zJy. Our limit to 
pulsed emission from CTB 109 2.3 mJy at 600 MHz corre- 
sponds to about 500 /xJy when scaled to 1489 MHz, again 
assuming spectral index of -1.6. Another example is 3C 
58, a Crab-like plerion containing an X-ray point source 
which is higly suggestive of a young pulsar (Helfand et 
al. 1995 and references therein). Our limit of 1.1 mJy to 
pulsed emission at 600 MHz corresponds, after scaling, to 
about a factor of two larger than the present best upper 
limit of 0.15 mJy at 1400 MHz (Frail & Moffet 1993). 

One of the supernova remnants observed during this 
survey, G78. 2+2.1, also known as 7-Cygni, contains 
the bright 7-ray source 2EG J2020+4026. Brazier et 
al. (1996) have recently discovered an X-ray source within 
the error box of 2EG J2020+4026. Assuming this to be the 
X-ray counterpart, they propose, on the basis of the 7- 
ray to X-ray flux ratio, that this source is most probably a 
Geminga-like neutron star. No significant radio pulsations 
were detected during our search of this remnant and our 
estimated 400 MHz lower luminosity limit of 1 1 mJy kpc 2 
certainly rules out the presence of a bright, favourably 
beamed, Crab-like radio pulsar associated with this source 
of emission. 



6. Pulsar Supernova Remnant Pair Statistics 

In this section, we address the question: How many pulsar- 
supernova remnant pairs (which we shall hereafter refer to 
simply as pairs) are likely to occur by chance in the present 
sample? As mentioned previously, this question has been 
tackled in some detail by Gaensler & Johnston (1995a&b) 
who concluded that the majority of claimed associations 
are likely to be chance alignments. The main motivation 
for the present analysis is partly to approach the question 
from a different direction than the modelling of Gaensler 
& Johnston (1995b). Our analysis has the main advantage 
that it is somewhat simpler and more model-free than the 
method described by Gaensler & Johnston (1995b). 

For the purposes of the analysis it is useful to charac- 
terise each pair by the dimcnsionlcss parameter (3 which is 
independent of distance, defined as the ratio of the angular 
separation between the pulsar and the remnant centroid 
to the angular radius of the remnant. Thus pairs in which 
the pulsar is within the boundary of the remnant occur if 
(3 < 1, whereas (3 > 1 indicates that the pulsar lies outside 
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Fig. 3. Top panel: The distribution of (3 for the sample of pulsar and supernova remnants after applying shifts in Galactic 
longitude to the pulsar sample; this is in excellent agreement with the distribution dN oc /3d (3 shown by the solid line. The 
lower panel shows the median characteristic age of the pulsars in the pairs as a function of /3 (see text). 
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Fig. 4. Top panel: The observed distribution of f3 for the sample of pulsar and supernova remnants. When compared to the 
distribution expected by chance shown by the solid line, this shows a clear excess of pairs with (3 < 1. The lower panel shows 
the median characteristic age of the pulsars in the pairs as a function of (3 (see text). 



the remnant (see also Shull et al. 1989; Frail et al. 1994; 
Gaensler & Johnston 1995 a&b). 

We are interested in deriving the distribution of [3 that 
occurs by chance and comparing this directly with the ob- 
served distribution. For completely unrelated sets of pul- 
sars and supernova remnants, the number of pairs occu- 
pying an annulus between /3 and f3 + d/3 is proportional 
to (3, regardless of the relative densities of pulsars and 
supernova remnants over the plane of the sky. To demon- 



strate this, we decoupled the respective pulsar and super- 
nova remnant samples by applying a systematic shift to 
the Galactic longitude of each pulsar and then calculat- 
ing the distribution of (3. To improve the statistics, we 
performed this procedure for shifts of ±4 and ±8 degrees 
in the Galactic longitudes of the pulsars. The shift sizes 
were chosen so that they are small compared to changes in 
the density of both types of objects on the sky, whilst be- 
ing larger than the angular size of any supernova remnant 
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in Green's (1996) catalogue. The shifted samples therefore 
contain only pairs which are truly unrelated, allowing us to 
deduce the expected distribution of chance remnants. This 
distribution is shown in Fig. 3 and is in excellent agree- 
ment with the theoretical prediction dN oc (3 d(3 shown by 
the straight line fit through the origin. The slope of the 
best-fit straight line shown in Fig. 3 is 47 ± 3. 

We are now in a position to apply this relationship, 
which is based on a superposition of four longitude-shifted 
samples, to the observed distribution of (3 shown in Fig. 4. 
The straight line shown in this case thus has a slope which 
is one quarter the value derived above. Comparing the 
observed distribution with this straight line we see a clear 
excess of pairs with (3 < 1 in the observed sample, whilst 
the distribution with (3 > 1 is in good agreement with the 
theoretical prediction. The pair excess is clearly due to 
the fact that many of these are genuine associations, not 
chance line-of-sight alignments. From Fig. 4, we infer that 
12 ± 5 pairs are likely to be genuinely associated. 

The difference between the shifted and observed sam- 
ples can also be seen in the lower panels of Figs. 3 and 4, 
where we have plotted the median characteristic age of the 
pulsars as a function of (3. The shifted samples show no 
significant deviation from a flat distribution with respect 
to (3, and have a median age of > 1 Myr. The observed 
sample however is clearly more youthful for (3 < 1 than 
for f3 > 1 for which the median characteristic age is again 
> 1 Myr. 

From an inspection of the observed sample and the 
available literature, we have compiled a list of the most 
likely associations in Table 4. A critical appraisal of each 
of these and other proposed associations can be found in 
the original references and see also Gaensler & Johnston 
(1995b) and Kaspi (1996). Note that this list is concerned 
only with Galactic associations and therefore does not in- 
clude PSR B0540— 69 associated with the supernova rem- 
nant in the LMC (Seward et al. 1984). Based on the above 
analysis and the available statistics, which suggested a to- 
tal of 12 ± 5 pairs, we expect anything from zero to nine 
further pairs to represent genuine associations. Examples 
of further candidates include recently discovered remnants 
around PSRs B1643-43 and B1706-44 (Frail et al. 1994), 
PSR B1930+22 and G57.3-1.2 (Routledge & Vaneldik 
1988) and PSR J0538+2817 and G180.0-1.7 (Anderson 
et al. 1996). In the absence of further information to vali- 
date these associations, we list the 8 sources in Table 4. 

As discussed in Sect. 5.1, both of the newly discovered 
pulsars in this survey, PSRs J0215+6218 and J1957+2831, 
have characteristic ages much larger than that expected 
for their target remnants. If they were to be genuinely as- 
sociated with these remnants, then both the pulsars would 
have had to have been born with periods similar to the 
presently observed values, 0.5 and 0.3 seconds, in order to 
explain the large characteristic ages. In addition, the re- 
spective observed value of (3 for these pulsars 0.7 and 0.8 
does not place them in the group of pairs most likely to 
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Table 4. A compilation of the most likely genuine pul- 
sar-supernova remnant pairs with /3 < 1. Our statistical anal- 
ysis suggests that, at most, a further nine pairs may be real 
(see text). 



be genuinely associated. Based on the results of the above 
analysis, there is no strong statistical requirement for ei- 
ther of these pulsars to be associated with the remnants 
although they cannot be entirely ruled out. 

Whether a significant number of radio pulsars were 
born with such long periods is controversial. Several au- 
thors, notably Vivekanand & Narayan (1981), Narayan 
(1987), Narayan & Ostriker (1990) and Deshpande et 
al. (1995) have found evidence for "injection" of pulsars 
into the population with periods ~ 500 ms, however 
other authors (Lyne et al. 1985; Stollman 1987; Lorimer 
et al. 1993) find no requirement for it. From an inspection 
of the pulse periods listed in Table 4, with the possible 
exception of PSR B2334+61, none of the pulsars listed in 
Table 4 are likely to have had such long periods at birth. 
Indeed Kulkarni et al. (1993) argue that PSR B2334+61 
is the energy source to G114.3+0.3, which appears to 
be a Crab-like nebula. In this case, the initial period of 
B2334+61 must have been < 100 ms. Thus the simplest, 
and most likely conclusion, to be drawn from the present 
sample of pulsar-supernova remnant pairs is that they 
support the notion that all pulsars are born with initial 
spin periods ^ 100 ms. 

Finally we wish to point out that, whilst our method is 
in principle sensitive to pairs over a large range in (3, the 
sample that we have used for our analysis is far from ho- 
mogeneous. Therefore any real pairs with (3 > \ are most 
likely to be underestimated in the present sample in com- 
parison with those with (3 < 1 since many of the targeted 
searches for pulsars in remnants have concentrated mainly 
close to the remnant centre. Indeed, together with the 
searches by Gorham et al. (1996) and Kaspi et al. (1996), 
our search represents the first major effort to search the 
entire area, rather than just the centroid, of many of the 
more extended supernova remnants. Further searches with 
improved sensitivity to pulsars with /3 > 1 will improve the 
situation. New multibeam searches of the Galactic plane, 
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presently underway at Parkes and Jodrell Bank, should 
provide a more homogeneous sample for further statisti- 
cal studies. 

7. Conclusions 

We have conducted a sensitive search for pulsars in su- 
pernova remnants. The search detected a total of three 
pulsars, two of which were previously unknown. The new 
pulsars, J0215+6218 and J1957+2831, were found dur- 
ing searches of the supernova remnants G132.7+1.3 and 
G65. 1+0.6 respectively. The case for associations between 
the new pulsars and their target remnants is presently 
unclear. From an analysis of the present sample of radio 
pulsars and supernova remnants, we reach two main con- 
clusions: (i) The number of real associations in the present 
sample is at most 12 ± 5. (ii) These are most likely to oc- 
cur for pairs with (3 < 1 for which we see a clear excess 
in the observed distribution compared with that expected 
by chance. 

Based on this analysis, neither of the two newly dis- 
covered pulsars seem likely to be genuinely associated. In 
order to unambiguously confirm/refute both these candi- 
date associations, measurements of the pulsar proper mo- 
tions are required. Such measurements would determine 
in each case whether the pulsar velocity is of the correct 
magnitude and direction to carry them to their presently 
position with respect to the remnant centre. Interfero- 
mctric measurements to determine the proper motion of 
both pulsars using MERLIN, the Multi-Element Radio 
Linked Interferometer network operated by the University 
of Manchester have recently been initiated. 
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